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Abstract 

There has been recent progress on computing real-time equilibrium 3-point functions in finite- 
temperature strongly-coupled AA=4 super Yang- Mills (SYM). In this paper, we show an example of 
how to carry out a similar analysis for a 4-point function. We look at the stopping of high-energy 
"jets" in such strongly-coupled plasmas and study the question of whether, on an event-by-event 
basis, each jet deposits its net charge over a narrow (~ l/T") or wide 1/7") spatial region. We 
relate this question to the calculation of a 4-point equilibrium correlator. 



I. INTRODUCTION 



Real-time, equilibrium, retarded Green functions are important to the study of relaxation 
in finite-temperature field theory. Through the fiuctuation-dissipation theorem, they are di- 
rectly related to the dynamics of how small deviations from equilibrium relax. In the context 
of strongly-coupled gauge-theory plasmas with gauge-gravity duality, 2-point retarded Green 
functions have been used to extract viscosity and other hydrodynamic transport coefficients 
[1, 2]. 3-point retarded Green functions have been used to (i) extract the stopping distance 
for certain types of high- momentum "jets" created in the plasma [3-5], and (ii) reproduce 
earlier results for second-order hydrodynamic coefficients [6, 7] and find the first correction 
to the large-coupling limits for such coefficients [7]. Finding useful analytic expressions for 
generic finite-temperature Green functions can be a great challenge, but problems of inter- 
est are often limiting cases where the Green functions simplify. For instance, hydrodynamic 
transport coefficients only depend on the low-momentum limit of the 2-point (or 3-point) 
function, whereas the 3-point function relevant for the jet problem is one where two of the 
three momenta are very high and the other is low [3]. In this paper, we will push the 
analysis of finite-temperature Green functions in such theories one point higher: we will 
discuss a problem, related to jet stopping, that requires the analysis of a 4-point function. 
In this application, the four momenta involved will all be either high or low compared to 
the temperature, which will make the analysis tractable. 

Here is the problem we study. We will loosely refer to spatially localized, high-momentum 
excitations of the plasma as "jets." Formally, one way to create such jets is to briefly turn 
on localized, high-momentum source terms in the field theory. In our earlier work [3],^ we 
tracked the progress of the resulting jet by measuring the subsequent evolution {j^{x)) of 
the density j° of a conserved charge carried by the jet. We found that the jet's charge 
was deposited along the jet's trajectory according to a distribution depicted in fig. 1. The 
charge subsequently diffuses hydrodynamically. In fig. 1, L is the space-time size of the 
localized source that we used to create the jet. Computing however, only gives us 

the averaged behavior of a jet. From this result alone, one does not know what happens 
on an event-by-event basis. Does each jet (i) deposit over a wide range of distances? Or, 
in contrast, does each jet (ii) stop and deposit its conserved charges in a small localized 
region, with fig. 1 representing a probability distribution for how far the jet travels? Later 
work [4], which provides a simple way to reproduce fig. 1 in terms of the trajectories of 
classical 5-dimensional particles in AdS5-Schwarzschild, supports the second interpretation. 
It should be possible to verify this interpretation by examining correlations {j'^ {x) {y)) of 
the charge density in the presence of the source. If charge deposition is indeed localized 
for each individual event, then there should be no significant correlation between charge 
deposition at significantly different distances. 

As we shall review, the evolution (i°(a;)) of the charge density in the presence of the 
source is related to an equilibrium, retarded 3-point function between (i) the measured 
charge density and (ii) source operators associated with the creation of the jet. The cor- 
responding Witten diagram in the gravity dual is depicted in fig. 2a. Correlations such as 
{j^ {x) {y)) will correspond to equilibrium 4-point functions, and fig. 2b gives an example 
of a corresponding Witten diagram.^ Our goal in this paper is to analyze such diagrams. 



^ See ref. [5] for a very short, breezy, and detail-free overview of ref. [3], 

^ For a recent discussion of 4-point current correlators at zero temperature, see ref. [8]. 
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FIG. 1: A qualitative picture of the average deposition of charge as a function of for jets created 
by the source described in Sec. II B and in ref. [3]. 

We will see later that, for the question of interest here, the internal line of fig. 2b is a bulk- 
to-bulk propagator with both bulk points located near the horizon. There is a relatively 
simple formula for the near-horizon bulk-to-bulk propagator, which will be a key feature in 
making the analysis tractable. 

At the end of our analysis, we will see that there is indeed no significant correlation 
between charge deposition at positions separated by ^ 1/T. 

In the next section, we discuss exactly what field theory correlator we need to measure 
in order to look for correlations between jet charge deposition at different locations. In 
particular, we will discuss how to define a correlator {Q{x) Q{y))jet of "charge deposition" 
Q{x) and what refinements are necessary to only capture the physics that we want. Section 
III is devoted to a discussion of exactly what type of equilibrium 4-point correlator we need in 
order to evaluate {Q{x) Q{y))jet and its relatives and the corresponding Witten diagrams. We 
start with a review of (r,a) notation for Schwinger-Keldysh formalism. Then we show that 
we need a 4-point correlator known as Gaarr? and we argue that there is only a single Witten 
diagram (fig. 2b with certain specifications for the causality properties of the propagators) 
relevant to the jet deposition question posed in this introduction. In preparation for the 
calculation of that diagram, section IV examines the equilibrium bulk-to-bulk correlator Grr 
in the gravity dual, which will be needed for the internal line in fig. 2b. We show how it 
can be related to standard bulk-to-boundary propagators. It turns out that only the case 
where both bulk points are very close to the horizon will be relevant, and we show that the 
propagators all have a simple structure in that limit. With all these preparations in hand, 
in section V we finally calculate the diagram and show that charge deposition is localized on 
an event-by-event basis. Various matters are left for appendices, and in particular Appendix 
B relates our discussion of bulk-to-bulk G^r to a recent discussion by Caron-Huot, Chesler, 
and Teaney [10]. 

II. GENERAL SETUP 
A. Notation 

We will use capital roman letters to denote 5-dimensional space-time coordinates (e.g. 
/, J = 0, 1, 2, 3, 5) and lower-case Greek letters for 4-dimensional space-time coordinates 
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FIG. 2: Examples of Witten diagrams in AdSs-Schwarzschild space for (a) a 3-point correlation 
contributing to {j^{x)) in the presence of source operators at xi and X2, and (b) a 4-point correlation 
similarly contributing to {j^ (x) (y)) . In both cases, the first diagram depicts traditional Witten 
diagrams drawn in imaginary time, with the boundary {u=0) represented as a circle. The second 
diagram depicts the same diagram in real time superposed atop a Penrose diagram of AdSs- 
Schwarzschild. For the type of correlators that we will examine in this paper, the bulk vertices in 
the latter case are confined by causality to the right-hand quadrant. Red lines indicate propagators 
that have large 4-momenta in our calculation (coming from the source) if one transforms to 4- 
momentum space, and blue lines (associated with measuring hydrodynamic response of charge 
densities) have low 4-momenta. We have used artistic license when drawing the boundary of the 
Penrose diagram with all four sides straight [9]. 



Throughout this paper, we generally refer to the fifth dimension of AdSs as u, where 
the boundary is at m = and the bulk has u > 0. Much of the discussion in the paper 
will be general but, when we get down to some specifics later on, we will choose u with 
AdSs-Schwarzschild metric 



(e.g. /i,z/ = 0,1,2,3). 




(27rT)2 



u 



{-f dt^ + dx^) + —du^ , 



(2.1) 



where R is the radius of AdS space, T is the temperature. 



(2.2) 
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and the horizon is at m = 1. At that time, we will also work in units where 2'kT — 1. 
For 4- vectors V^, we define light-cone coordinates by 

V^ = V^±V^, 14 = iyT ^ i(y3^y0)_ ^2.3) 
When writing integrals over 4-momenta, wc will use the shorthand notation 

B. The Source 

One could study this problem by creating a "jet" with most any localized high-energy 

source and then tracking any conserved charge density. The calculation is easier with some 
choices than others, and here wc will follow ref. [3] and study the propagation of R charge 
density, using a source that creates R charge. 

To later motivate our calculation of jet charge-charge correlations, it will be convenient to 
discuss sources that create jets with R charge ^ 1. On the other hand, it is also convenient 
to keep the detailed structure of the calculations here as close as possible to our previous 
calculations in ref. [3], where the source created jets carrying a single unit of R charge. As 
a compromise, in the main text we will use the single-charge source of ref. [3], described 
below, when writing down specific formulas involving the source. But little depends on this, 
and in Appendix A, wc discuss how the formulas generalize to cases where the R charge of 
the jet is large compared to 1. 

So, following ref. [3], we modify the 4-dimensional field theory Lagrangian by 

+ (2.5a) 

where are the SU(4) R-charge currents of the theory and Ac\ is a classical external source. 

Wc choose the external source to have the form of (i) a high-energy plane wave e^^'^ times 
(ii) a smooth, slowly varying, real- valued envelope function Kl^x) localizing the source to a 
space-time region of size L. Exactly as in ref. [3], we take 

= s^Ma e^'-" + h.c] A^(x), (2.5b) 

where 

F = (E,0,0,E) (2.6) 

is a very large light-like 4-momentum with frequency E ^ T\ Ma is an arbitrarily small 
source amplitude; e is a transverse linear polarization, such as 

£'^ = (0,1,0,0); (2.7) 
and are Pauh matrices for any SU(2) subgroup of the SU(4) R-symmetry, with 

r^ = T^±ir2. (2.8) 
A simple example of an appropriate envelope function would be 

Ai(x) = e-'^{^o/Lf^-\{x^/L)\ ^2.9) 
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L should be chosen large compared to 1/E, so that the momentum components in the source 
are all close to (2.6), but small compared to the large stopping distance that we wish to study. 
The creation of R charge by the source (2.5) is analogous to a decay ud of a. very 

high momentum W boson creating a localized excitation with isospin in a standard-model 
quark-gluon plasma. 



C. The Measurement: Previous Work 

In previous work [3], we measured the late-time behavior of 

{j^nx))A^, (2.10) 

for a system that started in thermal equilibrium. The superscript "(3)" indicates the R 
current associated with t^/2 in the SU(2) subgroup referenced by (2.8). The subscript "A^" 
indicates that the expectation is taken with the source term (2.5a) present in the Lagrangian. 
As reviewed in ref. [3], expanding to leading order in the small-amplitude source gives 

{j^'H^))a^, = i J d'x, d'x, Gi"^3)"^'^(xi, X,- x) Al,,{x,) 4,,i(a:2), (2.11) 
or equivalently 

= \ [ G^^^'^'^^^iQu Q2; Q) AZ,{Qi) 4:01(^2) e^«-(27r)^(5(^)(Qi + Q2 + Q), 
-'<3iQ2Q 

(2.12) 

where Gaar is the retarded, equilibrium 3-point Green function associated with source oper- 
ators j ■ A^i at Xi and X2 and the measurement operator j(^)'' at x. 

In this paper, we will use (r,a) notation to indicate different orderings of operators in 
n-point thermal Green functions. We review and summarize this notation (and fix our nor- 
malization conventions) in section III A, but for the moment it's enough to remark that 
G^aa -ar(2^i, 3^2, • • • , Xn-i] x) denotes an n-point retarded Green function where x is the mea- 
surement point and Xi, ■ ■ ■ , are source points. 

After a high energy jet stops and thermalizes in the plasma, its charge density diffuses 
out from the place where the jet stopped. What we most directly want to know is where the 
charge density was deposited [3, 11]. If locally thermalized charge is deposited with density 
Q{x) in the plasma, the late time evolution of the charge will be given by the diffusion 
equation 

{dt-DV^)f{x) = e{x), (2.13) 

where [12] 

is the R-charge diffusion constant. We can turn (2.13) around into an operator definition of 
the charge deposition:^ 

Q{x) = {dt-DV^)f{x). (2.15) 



The QQ{x) of ref. [3] is {Q{x))a^i in our notation here. 
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(b) 



FIG. 3: A cartoon of the space-time distribution of (a) {Q{x)) and (b) {T,q(x)) in the region of 
where there is net charge deposition. (See the discussion of fig. 9 of ref. [3] for a more detailed 
description.) 



Here and throughout this paper, we will often denote the measurement operator j(^)'* by 
simply when there is small chance of confusion. 

The charge deposition function G)(x) is only significant very close to the line a:° = of 
travel of the high-energy jet. We shall not concern ourselves with the detailed profile of the 
width of this function. And so, following ref. [3], we will write 

{e{x))A^, ^ h{x-) (Se(t=oo, x))a,„ (2.16a) 

where ^ 

Ee(x) = f dt'e{t',x), (2.16b) 

J — oo 

and the L subscript on 5l{x~) indicates that the 5-function approximates a function whose 
width is of order L. 

See ref. [3] for qualitative discussion of what the functions ©(x) and Se(^) look like. 
Here we will just summarize that the fact that net charge deposition is restricted to a region 
near the light cone, as shown by the cartoon picture of {Q{x)) in Fig. 3a, means that the 
time-integrated charge deposition (IIqIx)) will be constant in time for times well inside the 
light-cone, as depicted in Fig. 3b. We can summarize this constancy as 

{tQ{t=x'+e, x'))a^, ^ {te{t=oo, x'))a^, (2.17) 

whenever e ^ L. The important thing to remember is that even though we will find it 

convenient to talk about and calculate TiQ{t=oo,x^), that quantity is determined by charge 
deposition that takes place at t ~ x'^, not at t = oo. 

The great advantage of phrasing the problem in terms of the t^oo limit of T,q[x) is that 
the calculation in the gravity dual turns out to then only depend on the behavior of the 
corresponding 5-dimensional excitation at late times, when it has fallen very close to the 
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event horizon of the black hole. For the purpose of calculation, 5- dimensional propagators in 
Witten diagrams turn out to be much simpler and easier to deal with when the bulk points 
are very close to the horizon. 

It's important to note that most of the time, a localized, small- amplitude source will 
not have any effect at all on the plasma. It is only on rare occasions, proportional to the 
square \ A\^ of the source amplitude, that the source will create an excitation. When it does, 
the excitation will have the quantum numbers of the source. The jets will have energy and 
momentum ~ E. Their R charge ejet will be the R charge of the source operator,^ which is 
ejet=l for the specific choice (2.5b) and may take other values for the examples discussed in 
Appendix A. The average charge created in response to the source. 



A,: 



(2.18) 



is therefore ejet times the probability Vjet that a jet is created. In general, we will be interested 
in average charge densities and charge deposition in the case that a jet has actually been 
created. We can factor out the probability of creating the jet in the first place by defining 
expectations 

(•••)m^^^-^7T77^- (2.19) 

/ jet Wi/ Cjct 

So, in the notation of this paper, the average charge deposition of our jets is given by 

(e(a;))jet = ^|^. (2.20) 

Finally, at a practical level, the calculation of the Green function in (2.11) took advan- 
tage of approximations based on the momentum scales relevant to the problem. Since we 
measure where the charge is deposited by its subsequent diffusion, and since diffusion is a 
hydrodynamic process that takes place on distance and time scales large compared to the 
mean- free path and so large compared to 1/T, we do not need to resolve the structure of 
(j(^)°(x)) on scales as small as 1/T. In particular, we do not need to know the integrand of 
(2.12) except for when the components of Q (the 4-momentum conjugate to the measure- 
ment position x) are small compared to T. In contrast, the momenta Qi and Q2 associated 
with the source factors in (2.12) are both large compared to T because we have chosen 
a high-energy/high- momentum source. That means that, after Fourier transforming the 
boundary points, the (red) propagators associated with the source points xi and X2 in the 
Witten diagram of fig. 2a will have high 4-momcntum, and so can be treated in a WKB-like 
approximation, whereas the (blue) propagator associated with the measurement point x 
will have low 4-momentum, which also allows for simplification. Similar simplifications will 
occur when calculating 4-point correlators such as fig. 2b. 



Here and throughout, when we talk about the R charge ejet of a jet, we mean the charge corresponding 
to the specific U(l) subgroup of SU(4) that we have chosen for our measurement operator, i.e. the charge 
corresponding to what we have labeled j(3)o (2.10) and (2.18). 
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FIG. 4: A qualitative picture of charge deposition in a single event if charge deposition is very 
localized on an event-by-event basis. In this case, the location I of where the charge is deposited 
would have a probability distribution given by fig. 1. 



D. The Measurement: This Paper 

1. A first attempt 

To motivate our calculation, imagine that we used a source that creates jets with large 
R charge (ejet ^ !)• Then ask whether, in a single typical event, the spatial distribution of 
deposited charge looks like (i) fig. 1 or (ii) fig. 4. In the latter case, fig. 1 would give the 
probability distribution for the position / in fig. 4. Naively, a simple way to distinguish these 
two cases would be to measure the average correlation 

(e(^)e(^)),.= ^ Q/e'! (2.21) 
for well separated x and y. The idea is that, for well separated x and y, case (i) would give 

(e(a;)e(?/)>.^,^(e(x)>^^X0(i/)>M' (2.22) 

whereas case (ii) would instead give 

(0(x)0(y)>^^^^O (2.23) 

if X and y were much further apart than the width the peak in fig. 4. 

When discussing correlators, there are operator ordering considerations except in the 
classical limit. Note, for instance, that the right-hand side of (2.22) is (i) symmetric with 
respect tox ■H- y and (ii) real, whereas the left-hand side is not precisely either. Factorization 
such as (2.22) would require the imaginary part of the correlator to be small. The real 
and imaginary parts of (0(x) 0(y)) are given respectively by the symmetric combination 
^{{Q{x), Q{y)}) and the commutator —i^{[Q{x), Q{y)]). In this paper, we will focus on the 
real part, for which the two cases distinguished above would be 

i({e(x),e(y)}>^^^^(0(x)>^^Xe(y)>^^^ (2.24) 

and 

i({e(a;),e(y)}>^^^^0. (2.25) 
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As we will sec in a moment, the idea for our correlation measurement will require a little 
bit of refinement, but the basic goal of this paper is to demonstrate a refined version of (2.25) 
and so verify that the deposition of charge by a jet looks like fig. 4 on an event-by-event 
basis with a width Ax^ <1/T. We will refer to this as "localized" deposition of charge. 

Eqs. (2.24) and (2.25) lay out two logical extremes for how a correlation might behave, 
but a priori the truth could instead lie somewhere in between. In particular, factorization 
(2.24) is not a very plausible option for non-local correlations unless ejet^l- But the simple 
question of whether the correlation is local or non-local can be sensibly asked for ejet=l 
as well as ejet^l- In the main text, for the sake of concreteness and compatibility with 
previous work, we will focus on the case ejet=l corresponding to the specific source (2.5b). 
For the sake of completeness, we will discuss the same question for ejet^l in Appendix A. 
In both cases, we will find locality for the suitably-refined correlation. 

In calculations, we will find it easier to study the correlator 

i({E0(x°=oo, x), S0(/=oo, y)}) (2.26) 
of time-integrated charge-deposition (2.16b) rather than ^{{Q{x),Q{y)}) directly. 

2. Refinement 

When we later evaluate the correlation |({0(a;), 0(y)}) in the gravity dual, we will see 
that we want to throw out certain types of diagrammatic contributions which do not have 
the causal structure appropriate for the question of how charge deposition from a jet is 
correlated. In this section, we discuss what goes wrong if one includes all contributions to 
|({G(x), ©(?/)}), and we outline a more careful definition of charge deposition which avoids 
unwanted contributions and only includes the physics of interest. 

From our definition (2.15) of charge deposition Q{x), the correlator we proposed above 
for study is 

l{{e{x),e{y)}).^^ = {d^o - DVl){dyO - DVl)\{{3\x),f{y)}).^^. (2.27) 

This would be fine if the only source of charge in the system contributing to this cor- 
relator were the charge deposited by the jet. However, even if there were no jet at all, 
the equilibrium correlator {{j^{x)^j^{y)})c(i would be non-vanishing and would measure the 
correlation of charge-density fiuctuations in the plasma. Not only is the equilibrium cor- 
relation not the physics we want to measure for studying jet stopping but, unlike the case 
of charge deposition, its time dependence is not the simple retarded time-dependence of 
the diffusion equation implicitly assumed in the introduction (2.13) of the charge deposi- 
tion distribution Q{x). To see this, consider that the equilibrium charge density correlation 
function {{j^{x),j^{y)})eci can be related to the retarded charge density correlator via the 
fluctuation-dissipation theorem:^ 

K{/'/}>eq = -cth(» ImG°°(a;,qr), (2.28) 



^ See, for example, the discussion by Kubo [13]. Eq. (2.28) above is a special case of our later (3.3). 
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where we work in 4- momentum space {u, q). In the hydro dynamic hmit of small frequency 
and momentum, G'^ oc Dq^ /[—iuj + Dq^) has retarded diffusive behavior, and (2.28) gives 



where x is the charge susceptibility.^ It's useful to also write this as 

Dq^ Dq^ 



-icu + Dq^ icu + Dq^ 



(2.30) 



There is not only a diffusive pole —iu + Dq^ — but also a conjugate pole iu + Dq^ — 

corresponding to its time reversal. Multiplying both sides by factors of —icu + Dq^ to turn 
j°'s into 0's would simplify the first term of (2.30) but not the second. The difficulty is that 
the correlation represented by {{j'^^j^})^^ is time symmetric, unlike the diffusion operator 
dt — -DV^ and unlike the physics of charge that is suddenly deposited in the medium by a 
thermalizing jet and only subsequently diffuses. The upshot is that the equilibrium contri- 
bution to ({i°(2^),i°(y)})eq messes up our desired physical interpretation of the correlation 
(2.27). 

There is a simple fix to the problem identified so far: Subtract out the equilibrium 
contribution to the charge correlation, and so study 

(9.0 - DVl){d,o - DVD [<0'°W> /(?/)»jet - ({/(^)>/(y)}>eq] (2-31) 
instead of (2.27). 

But there is still a problem. Consider the case where the jet passes through a region of 
plasma that contains some equilibrium fluctuation in the charge density. It's conceivable 

that interaction with that fluctuation, depending on its sign, could bias the probability that 
the jet stops earlier or later, as depicted in fig. 5. So we could get a contribution to (2.27) 
where the j^{x), for example, represents a pre-existing charge fluctuation in the plasma (i.e. 
not something caused by the jet) and the j°(?/) represents diffusion from a jet depositing its 
charge at y. Let's call this a plasma-jet contribution to the correlation. Such a contribution 
would contaminate what we really want to know, which is the correlation between the 
jet depositing its charge at x and depositing its charge at y. Plasma-jet contributions to 
the correlation would also not give the retarded time dependence in a;° assumed by the 
{d^o — D'V^) factor in (2.27) for the same physical reason as in the discussion of the purely 
equilibrium correlator {{j^{x),j^{y)})cq. This vague concern will be made concrete when we 
later discuss diagrammatic contributions to (2.27) in the (r,a) formalism. We will sec that 
some diagrams correspond to the physics of the charge deposition correlation that we want, 
while other diagrams correspond to plasma-jet correlations that we are not interested in for 
the purposes of answering our question about how charge is deposited. 



^ The overall normalization of the right-hand side of (2.29) is determined by integrating both sides du>/2TT 
over all w and then taking q — > 0. After this procedure, the left-hand side becomes the equal-time expec- 
tation ((total charge)^)/volume, which defines the charge susceptibility %. For a more specific formula for 
the hydrodynamic limit of G^{lo, q) for jV^=4 SYM in particular, see, for example, eq. (79) in the review 
article [1]. 
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FIG. 5: A cartoon depicting a possible correlation between a charge fluctuation in the plasma at 
X and the deposition of the jet's charge at y. 



Since neither the original correlator |({G)(x), 6(y)}) nor the subtracted version (2.31) 
exactly capture the physics we want, how can we define something that does? The problem 
is that our definition of charge deposition as given by (2.15) for Q{x) is not really a good 
definition of charge deposition when we want to study correlators. We should define more 
carefully what it means for the jet to deposit charge somewhere. Here is one way to do 
that. Until the jet deposits its charge, we know from previous work that it stays close to the 
light-cone x'^ = x^, within a distance x~ ~ L (assuming L > 1/T). We wish to define the 
net charge density deposited by the jet in the vicinity of some point x~^ along the light-cone. 
Consider two 3-dimensional hyper-planar strips defined by 



X 



±5, 



plane ; 



\x — X 



< 



-plane; 



(2.32) 



corresponding to the dashed lines in fig. 6a. Choose Spiane and £piane so that L <^ 5piane <^ 
^piane ^ Ax^, wherc Ax~^ is whatever scale at which one wants to resolve differences be- 
tween stopping distances. Aa;^ should definitely be chosen small compared to the typical 
stopping distance {ELY^^ /T. Consider the space-time region bounded by the two planes 
and approximately define the charge deposited at x~^ as the difference between the charge 
leaving that region through the top plane and entering it through the bottom plane: 



Se(oo, x^=|x^ 



-plane 
1 



t- plane 



' 2 ■^planc 



X ' 



'^plane _L X^ -\- 5j 



plane " 



_ / X -\- (5plane 



X ' 



^planc " 



(2.33) 



where the dS integral in the first line is over both hyper-plane segments, with dSfj_ pointing 
outward from the region between, and V± is the area of the transverse dimensions x±7 



^ Remember that we have chosen to set up a transverse translationally invariant problem, so V± is just an 



12 




(a) (b) 

FIG. 6: (a) The hyper-plane segments (dashed lines) used to define Se corresponding to charge 
deposition near the light-cone points x and y. (b) An enlarged picture of one set of these plains 
showing how a charge fluctuation from the plasma might sneak between a pair of planes through 
one of the open ends and so be wrongly counted in the charge deposition measurement. 



The jet sneaks into the region, without crossing either plane, through the lower-left opening 
between the planes in fig. 6a. If it does not thermalize there, then it sneaks out again, 
uncounted, through the upper-right opening. If some or all of its charge thermalizes in 
the region, that deposition of charge will be counted by (2.33). This definition solves the 
problem of jet-plasma correlations because a pre-existing charge fluctuation in the vicinity 
of X will give a canceling contribution to (2.33) between passing through the lower and upper 
planes, as depicted in fig. 6a. There will be errors in this procedure because of edge effects: 
some fluctuation in the equilibrium charge density might sneak in one side and be counted, 
as shown in fig. 6b. But such effects will scale away as one makes ipiane larger and so can be 
disentangled. 

We use the bar over Ee in (2.33) to distinguish our refined version of T,q{oo,x^) from 
our original definition (2.16b). The bar can also crudely be thought of indicating that 
this definition averages Se(oo,x^) over a distance of order £piane, and so introduces some 
smearing. Such smearing is inconsequential as long as one keeps ^piane small compared to the 
differences in stopping distances that one wishes to resolve. Finally, we remind the reader 
that the quantity Se(oo, x^) arises from charge that is deposited near the light-cone x^ ~ x^ 
and not at x°=oo. See (2.17). 

Consider now yet another possible source of {{j^ (x) , (y)}) correlation. When the jet 
stops, it deposits energy and momentum into the medium, which can then propagate away in 
the form of sound waves. These sound waves locally increase and decrease the temperature 
and pressure as they pass by. A small change ST in local temperature will give a small, 



infinite (and uninteresting) normalization factor. Eq. (2.33) really defines the aj^-averaged time-integrated 
charge deposition function. We could have equally well defined the symbol Se throughout this paper to 
be the aj^-integrated (rather than averaged) time-integrated charge density, in which case there would be 
no 1/V± factors in (2.33). 
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FIG. 7: A jet travels along the light cone and then deposits its charge, which diffuses, as well as 
sound energy, which propagates. The local oscillation of temperature associated with the sound 
wave can then modify the local correlations between j^{x) and j^{y) at, for example, the points 
shown. Such a change in the correlation would contribute to (2.31) but is invisible to the refined 
measurement (2.33), which is defined specifically for x and y on the light cone and only involves 
at points close to the light cone, such as at the dashed lines. 



local change to the contribution to {{j^{x), j^{y)}) in that region from finite-temperature 
fluctuations of the plasma there, leading, for example, to a contribution 



(d^o - DVDidyO - DVl) {{Ax), Ay)}) 



' eq,T+<5T 



eq,r 



(2.34) 



to the subtracted correlator (2.31). We will later see a Witten diagram in the gravity dual 
calculation which naturally produce this sort of effect. Fortunately, our refined measurement 
(2.33) would eliminate this contribution since it only looks at charge density near the light 
cone and so is not affected by effects on the charge density far away, such as shown in fig. 7. 

The definition (2.33) is an example of how to define charge deposition, but it would be 
awkward and difficult to compute |({S0(cxd, x^), Se(oo, y^)) using this detailed formula in 
practice. However, we will see that its effect is to simply eliminate certain types of (gravity 
dual) diagrammatic contributions to the charge-charge correlation and to reproduce our 
naive |({S0(cxd, x'^), S(cxd, y^)}) for the diagrammatic contributions which survive. The 
upshot is that in practice we will just compute 



i({Se(oo,x),Se(oo,?/)}) 



jet 



relevant diagrams 



(2.35) 



We will see which are the relevant diagrammatic contributions in sec. Ill D. 



III. 4-POINT FUNCTIONS AND DIAGRAMMATIC SETUP 

In this section, we see precisely which type of 4-point equilibrium correlator we need to 
calculate and how it corresponds to an appropriately causal real-time Witten diagram in 
the gravity dual. But first, it will be convenient to review the (r,a) notation for organizing 
the classification of real-time correlators, based on Schwinger-Keldysh formalism. 
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FIG. 8: The Schwinger-Keldysh contour in the complex time plane. 
A. Field theory review of (r,a) notation 

(r,a) notation is nicely summarized in ref. [14], and examples of its application to 
calculations with gauge-gravity duality may be found in refs. [10, 15, 16]. One begins 
with Schwinger-Keldysh formalism for a system in thermal equilibrium, where the finite- 
temperature imaginary-time integration contour from to —iP is deformed into real time as 
shown in fig. 8.*^ Operators may be placed on either the upper or lower horizontal sections of 
the path, with operators labeled type "1" or "2" correspondingly, and the n-point functions 
are then path ordered. For example, consider the ra-point correlator of a single field 0{x). 
Then 

^'^G'i2212(a;i, X2, Xg, X^, X5) = (Tp[Oi(Xi) 02(0:2) 02{X:i) Oi(x4) 02{X5)]) 

= (r[o(x2) o(x3) o(x5)] r[o(xi) oix,)]), (3.1) 

where 7p is Schwinger-Keldysh path ordering and T and T represent ordinary time ordering 
and reverse time ordering respectively. For 2-point functions, Gu is the Feynman propagator 
Gp] G22 is related to Gp] and G12 and G21 are the Wightman correlators G< and 
The (r,a) notation consists of switching the basis of the labels 1 and 2 to^ 

0^ = 0i-02, 0r = l{0i + 02). (3.2a) 

We will define the corresponding n-point Green functions by 

i''-'G^,...^„{x,, ■ ■ ■ ,a;„) = 2^^-' {rp[0^,{x,) ■ ■ ■ O^M]), (3.2b) 

where each is either "r" or "a" and riy- is the number of "r" indices. Different authors 
use different factor of 2 conventions, and ours are those of Wang and Heinz [14]. For 2- 
point functions, Ga_^{0,x) is the retarded Green function G^^ix); Gi.a(0,a;) is the advanced 



^ People sometimes find it useful to deform the contour so that there is a vertical segment that drops from 
+too to +too — icr on the right-hand side of fig. 8, and then the vertical segment on the left-hand side 
drops from —too — io' to —too — if^- For instance, ref. [18] uses a — /3/2. But here we will stick with the 
traditional a = 0. 

^ In the literature on the Schwinger-Keldysh formalism, and Or are sometimes referred to, up to factor 
of 2 normalization conventions, as the "quantum" and "classical" fields Oq and Od respectively. See, for 
example, ref. [15, 19]. 
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Green function Ga{x); iG„{0,x) is the symmetric correlator ({0(0), 0(x)}) = iGy +iG<; 
and Gaa = 0. An equilibrium relation that we will find useful later on is the fiuctuation- 
dissipation theorem, 

iGrriou, q) = -2 cth{\^u) Im G^{u, q) = cih{\^u) [iG'R(a;, q) - iGA{uj, q)] (3.3) 

for bosonic correlators, ^° which we have previously referenced in (2.28). 
Here is a useful property of correlators in (r,a) notation: 

Rule 1. An n-point function vanishes if the largest time max(ii, • • • , t„) is asso- 
ciated with an operator 0^. that has an "a" label. 

This follows because, if x° is the largest time in the correlation, then there is no difference 
between the operator ordering of 0i(x) and (p2{x), and so (f)3,{x) defined by (3.2a) vanishes. 
Rule 1 has the following corollary: 

Rule 2. The n-point function Gaa-.-a with all "a" labels is zero.^^ 

It is also useful to take note of a property of perturbation theory in (r,a) formalism, even 
though the 4-dimensional field theory we are interested in (strongly-coupled Af—A super 
Yang-Mills) is not weakly coupled. In the Schwinger-Keldysh formalism, using the time 
contour of fig. 8, the path integral schematically takes the form 

J [P0i(a;)][P02(a;)]e'^-(^['^il-^['^^l)e-^^['^3]^ (34) 

where 0i and (1)2 live on the upper and lower horizontal contour and 03 lives on the vertical 
stub of the contour at Ret= — too- In Hamiltonian language, the 0i and 02 parts respectively 
generate the evolution operators e'^S^'^^ and e'^^S^'^^ that appear in the time evolution 
^+ifHdtQ^-tfHdt operators that we wish to measure, while the 03 part generates the 
equilibrium density matrix e~^^ . Order by order in perturbation theory, there are generically 
corrections both to the evolution operators e^^^^*^^ and to the initial density matrix e~^^ . 
For our application, only corrections to the evolution will turn out to be important, and so 
focus on the corresponding interaction terms 

f (y[0i] - V[h]) = / {V[<l>r + i0a] - V[<l>, - |0a]). (3.5) 

J X J X 

The integrand is odd in 0a. As a result we have the following rule in (r,a) notation: 

Rule 3. Any interaction vertex (associated with time evolution) must have an 
odd number of "a" fields. 



In addition to Kubo [13], sec, for example, eq. (32) of Wang and Heinz [14], Also iG^{uj,q) = — 2[n(cj) + 
1] ImGR(w,q) and iG<{uj,q) = -2n(w) ImGR(w,q), where n{Lo) = iKeP"^ - 1) = i[ctli(i/3a;) - 1] 
is the Bose distribution function. For a textbook reference, see, for example, eq. (7.6.41) of ref. [20], 
where their 2%" is the commutator ([0(x), 0(a;')]) = iG-^ — iGk and their S is the Wightman function 
i(j> = {0{x) 0{x')). By taking w — >■ —uj one gets the corresponding relation for and then adding G> 
and G< gives Grr- As far as the case uj=Q is concerned, we have implicitly assumed (O) = 0, and so O 
should be replaced by O — (O) if not. 

This rule is equivalent to the largest-time diagrammar identity of Veltman [21]. 
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FIG. 9: The non- vanishing (r,a) assignments for the propagators in one of the diagrams that 

contribute to Gaarr{xi,X2,X,y). 

As a simple example, consider the contribution of the diagram of fig. 9 to the 4-point 
correlator Gaa.rr{xi,X2,x,y). Applying Rules 1-3, the figure shows all the non-vanishing 
(r,a) assignments for this diagram's contribution to Gaarr- 

Of course, we will not be doing perturbation theory in 4-dimensional field theory for 
strongly-coupled J\f=4: super Yang-Mills theory. However, as others before us [10, 16], we will 
extend the use of the (r,a) notation to organize calculations in the 5-dimensional gravitational 
dual theory. In the gravitational theory, we are doing perturbation theory. Because the 5- 
dimensional AdSs-Schwarzschild background is periodic in imaginary time with period 27rT, 
the methods of thermal field theory can be taken over to this 5-dimensional field theory 
problem [17]. 

Finally, returning for a moment to ordinary perturbative thermal quantum field theory, 
we should mention that there is a long, involved story about whether and how, in the limit 
too oo, one may ignore vertices which correspond to interaction terms such as {(ps)^ that 
live on the vertical stub of the contour in fig. 8. (See Gelis [22].) However, this is a subtlety 
which may be blissfully ignored in many applications of interest. For example, consider a 
diagram like fig. 10, where the "3" vertex represents a (^s)^ interaction along the vertical 
stub of the contour in fig. 8. This diagram is zero because, by the same logic that led to 
Rule 1, the propagator G^a = 0. At leading order, there are no non- vanishing diagrams 
for Gaarr that iuvolvc corrections to the initial density matrix, and so all the vertices in 
fig. 9 just involve r's and a's. In the gravity dual, an analogous concern would arise if 
evaluating a correlator whose causality structure allowed bulk vertices to be inside the past 
or future horizon in fig. 2 and potentially touch the singularity, but this will not happen 
in our application. (When it does happen, ref. [23] shows how to evaluate correlators in 
a way that satisfies thermal equilibrium relations and does not require consideration of 
points inside the horizon. Alternatively, one may use equilibrium relations [14] to relate the 
correlator to another that avoids the problem. Or, as a matter of principle, one might use 
the complexified AdSs-Schwarzschild of Skenderis and van Rees [24, 25].) 

B. The desired 4-point correlator 

In our earlier work [3], we needed an equihbrium 3-point retarded correlator G'aar(a^i, ^2] x) 
to calculate the response (j°(a;))A<,i of charge density to the creation of a jet, where each 
of the "a" subscripts in Gaar is associated with the source operator. In this section, we 
will see that we need the equilibrium 4-point correlator Ga,a.rr{xi,X2;x,y) to calculate the 
correlator |({j°(a;), ^^(y)})^^!. The appearance of "r" subscripts for both x and y in Gaarr 
can be thought of as a generalization of the equilibrium relation iG„ = {{0{x), 0{y)})cq to 
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FIG. 10: An example of a diagram involving the fields (p^ on the vertical stub of the Schwinger- 
Keldysh time path of fig. 8. Here the vertex marked with 3's represent a {(ps)^ vertex on that 
stub. This assignment (and all other assignments of this diagram that involve (p^ vertices) vanishes 
because G^a = 0- 



the case of jets created by a source A^i- 

Write the Hamiltonian of our system as H{t) = Hq + 6H{t), where 6H{t) contains the 
small-amphtude source terms that we introduce to create the jet and Hq is everything else 
in the full, interacting Hamiltonian of the theory. Take 6H{t) to be localized in time, as 
previously discussed, and assume the system starts in equilibrium at t — )■ — oo. As reviewed 
in our earlier work [3], the response of an observable O is given by the following expansion 
in 6H: 

A{Oit)) = {0{t))H-{0)Ho = J rftiGar(ti;t) + ^ J dhdhG,Uti,t2;t) + -- - , (3.6) 

where the various Ga - ar are the equilibrium n-point retarded correlation functions, given in 
this case by 

zG,,{t,-t) = 9{t - h) {[Oit),6H{t,)])^^, (3.7) 
^'G,Uh, t2; t) = 9{t - t^) ^(ts - h) {[[0{t), 5H{t^)l5H{h)]) 

+ e(t - ti) e{U - h) {[\0{t), 5H{h)i 5H{t2)])^^, (3.8) 

etc. Taking O to be j^^^^(x), and SH from (2.5) to be / (Px j ■ Ad, reproduces the formula 
(2.11) that we reviewed earlier for {j^^^'^(x))A^i- 

Nothing changes in the derivation of the above formulas if one replaces 0{t) by a product 
of operators 0{t) 0{t') in the particular case that t and t' are both after the time at which 
the source has turned off (so that the source does not affect the time evolution operator 
between t and t'). Such late-time observables are precisely what we're interesting in for 
studying jet stopping. We will focus in particular on the symmetric product ^{0(t), 0(t')} 
of measurement operators. In the case of interest to jet quenching — a high-momentum source 
SH and low-momentum observables O — the leading contribution at post-source times to the 
6H expansion will then be 

lA{{0{t),0{t')}) ~ dhdt2G{t,,t2;t,t') (3.9) 



with 



t'G{h,t2; t, t') = e{h - ti) {[[l{0{t),Oit')}, SHit2)],5Hih)])^^ 

+ ih^t2). (3.10) 
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(a) (b) (c) 

FIG. 11: Witten diagrams contributing to the 4-point current correlator in our problem, with the 
sources at the bottom and the measurement points x and y at the top. For (a), one should also 
add the diagram where x and y are interchanged. The types of lines are described in the text. 

We have dropped factors such as 0{t — ti) and 6{t' — ti) because the formula (3.9) only applies 
in the case of post-source times, where t and t' are each greater than both ti and t2- In 
this case, it's easy to check that the (r,a) ordering rules described in section III A imply 
that G above is the same thing as iGg^^jl. The factor of ijl arises from our normalization 
convention (3.2) for (r,a) correlators. 

Taking the same choices for O and bH as before, we get the following 4-point generaliza- 
tion of the 3-point relation (2.11): 

iA({j(3)^(a;),j(^)^(y)}>^^^ = ^| d'x,d'x,G^:^f^-^^^{x,,x,-x,y)Al^J,x,)A\^^^^ (3.11) 

(for times after the source is turned off). 

C. Witten diagrams in the gravity dual 

Recall that we take our high-momentum source (2.5) to be translation invariant in the 
transverse directions (x^, x^) and transversely polarized. Restrict attention to measurements 
of correlators of the charge density i.e. /i=z/=0 in (3.11). The Witten diagrams which 
contribute to the desired 4-point correlator in (3.11) are shown in fig. 11. (Here, for the sake 
of pictorial simplicity, we have not bothered to draw Penrose-diagram versions as in fig. 2.) 
We are also not yet implementing the "aarr" prescription on Gaarr, which we will discuss 
shortly. In these diagrams, straight solid lines represent propagators associated with the 
transverse polarizations of the 5-dimensional gauge field. These have high 4-momentum 
and have index ± in the SU(2) subgroup of the SU(4) R symmetry, corresponding to (2.5b). 
Wavy lines represent the /=0 and 3 components of the gauge field Ai (in ^5=0 gauge), 
which mix and which we've abbreviated as A\^ in the figure. These wavy lines have low 
4-momentum (since they are associated with our late-time measurement of diffusion) and 
index 3 in the SU(2) R-symmetry subgroup. The dashed line can represent either a graviton 
or a massive scalar discussed by Romans [26]. 

As seen from fig. 11c, diagrams for the 4-point function can contain other supergravity 
(SUGRA) fields besides the ones dual to the boundary operators. The full set of fields in the 
SUGRA theory is rather complicated, but the problem of enumerating the diagrams can be 
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Gr(9) = " = 

q >• q »• 

Ga(9) = " = — ^— 

q >• q >• 

Grr(9) = ' ' = — H— 

9 *- 1 >• 

FIG. 12: Our diagrammatic notation for propagators in the (r,a) formalism. The same convention 
is used for both bulk-to-bulk propagators G^^^ and bulk-to-boundary propagators Q. 
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(a-2) 


(a-3) 


X 

(b-1) 


X 

(b-2) 
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V 




(c-1) 


(c-2) 





FIG. 13: Non-vanishing (r,a) assignments for Witten diagrams contributing to Gaarr in our problem. 
The first row is simply a rewriting of fig. 9 using the notation of fig. 12. 

simplified due to the fact that our sources all lie in an SU(2) subgroup of the full SU(4) R 
symmetry. Lii, Pope, and Tran [27] showed that the 5-dimensional SU(4) gauge theory in the 
gravity dual can be consistently truncated to a theory of just a SU(2)xU(l) subgroup, and 
the Lagrangian for A^=4 SUGRA in five dimensions with a gauged SU(2)xU(l) symmetry 
group is conveniently given by Romans [26]. However, none of the details matter because we 
will see later that only the diagram of fig. 11a will contribute to the physics we are interested 
in. 

Now focus on Gaarr in particular. In order to depict (r,a) assignments graphically, we will 
use the notation defined in fig. 12. Using the rules of section III A, the (r,a) assignments for 
the diagrams of fig. 11 which contribute to Gaarr are shown in fig. 13. For comparison, the 
diagram corresponding to the 3-point function Gaar is shown in fig. 14. 
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FIG. 14: Non-vanishing (r,a) assignment for the 3-point function G-, 



D. The dominance of diagram (a-1) 

We win see now see that only diagram (a-1) potentially contributes to the physics of 
interest to us. The other diagrams, which all involve a Q„ on one of the two measurement legs 
(the wavy lines), only contribute to plasma-jet and other correlations of the sort discussed 
and discarded in section II D 2. In all of the diagrams, we will make a low-momentum 
approximation for the 4-momenta Q and Q' conjugate to the measurement points x and 
y, as discussed at the end of section II C. That is, we will blur our eyes and only concern 
ourselves with resolving measurements on distance and time scales larger then 1/T. 

We start by discussing some of the properties of diagram (a-1). The small 4-momentum 
limit of the gauge bulk-to-boundary propagator Qf^ is given in refs. [1, 3].^^ The details do 
not matter yet except that the general form is 

where each fi^{u; u, q) is polynomial in oj and q (and so is local in 4-position space). Above, 
the pole {iuj — Dq^)^^ describes diffusion according to the diffusion equation. If we now 
construct O from j° for the measurement point, by multiplying (3.12) by iu — Dq^ as in the 
definition (2.15) of 9, we get 

{dt - DV^)gf^ = completely local in x. (3.13) 

That is, in 4-position space, {dt — D'V^)Qf^ connecting a bulk point {x',u) to a boundary 
point X only has support at x=x' (up to the 1/T smearing of resolution intrinsic to our small- 
momentum approximation for the measurements). So, once we apply the diffusion operators 
dt—D'V^ to both measurement points as in (2.27), diagram (a-1) becomes a truncated 
diagram depicted by fig. 15a, with some derivatives and possibly w-dependent factors that 
correspond to fi^{u] u, q) above acting on the wavy stubs. If we drew a corresponding 
diagram for our earlier calculation [3] of {Q{x))a^i, it would be fig. 15b. In that work, 
we found that the charge deposition only had support for x very close to the light cone, 
where close means within max(L,T~^). We will see the same thing happen in this paper in 
the calculation corresponding to fig. 15a. The underlying reason is that the 5-dimensional 
wave created by the source, given by a bulk-to-boundary propagator from a source point Xi 
(represented by a straight line from the boundary in figs. 13-15) convolved with the source 
wave packet flL{xi)e'''''^^ on the boundary, stays close to the x^=x^ 4-space light cone for all 
times relevant to charge deposition (x° < E^^^). We will leave the details of calculating the 



See in particular eq. (4.7) of rcf. [3 
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(a) (b) 



FIG. 15: Versions of fig. 13(a-l) and fig. 14 witli the measurement bulk-to-boundary propagators 
truncated by applications of the diffusion operator dt—D'V^. We have included a depiction of the 
boundary in these figures to emphasize that the truncated vertices lie in the bulk. 



contribution from diagram (a-1) until later, but this qualitative point about the truncated 
diagram will be enough to understand what happens with the other diagrams of fig. 13. 

In diagram (a-1) of fig. 13, the charge-charge correlator j^{y)}) a^i was given by 

deposition near the light cone, represented by fig. 15a, propagated forward in time to the 
measurement points x and y by diffusive propagators, contained in the directed wavy fines 
in fig. 13(a-l). We depict this situation pictorially in fig. 16(a-l). What happens if instead 
one of those bulk-to-boundary propagators for x or ?/ is a low-momentum Q„ instead of a 
Qr, as in diagram (a-2)? Applying (3.3) to the gravity calculation, 

Qrviu; q) = cth(|/3a;) [^r(m; u, q) - ^a(m; q)] 

= ci\i{\(3uj) [Gn{u; OJ, q) - g^iu; -u, -q)] . (3.14) 

As a result, ^rr('w;g) above is even under inversion q — t- —q of the 4-momentum. Fourier 
transforming in 4-space, that means that Q„{x' ,u] x) is even in the difference x — x' oi the 
boundary and bulk 4-positions. We depict this situation pictorially in fig. 16 (a-2), where the 
point associated with Q^r has a backward-time diffusive region, corresponding to the second 
term in (3.14), that is PT-symmetric with the forward-time diffusive region. As a result, if 
we calculate the refined time-integrated charge-deposition operator Se of (2.33) at a point 
on the light cone, the contributions from subtracting the integrals of j° over the two dashed 
hyper-plane segments in the figure will cancel. And so diagram (a-2) will not contribute to 
our refined measure ^{{Eq^oo, x^), T,q{oo, y^)}) of charge deposition correlations. The same 
conclusion holds for diagram (a-3) and, by similar logic, for (b-1) and (b-2). 
In what follows, we will abbreviate Se(oo,x^) as S(x^). 

Now we come to the (symmetrically related) diagrams (c-1) and (c-2) of fig. 13. Focus 
on (c-1) for the sake of specificity and let (x', u) be the location of the top-most bulk vertex, 
as indicated in fig. 17. Because of the fiow of causality (indicated by the arrows) from the 
source points Xi and X2, the bulk point {x',u) must lie in the 5-dimensional causal future 
of Xi and X2, and the measurement point y must lie in the future of that. This requires x' 
to be in the 4-dimensional causal future of the source and y to be in the causal future of 
x'. Now consider two cases: x' is (i) far from or (ii) near the x^ ~ x^ lightcone, as in figs. 
18a and b respectively. In the first case, it is causally impossible for the measurement point 
y to be near the light cone, and so there can be no contribution to the refined correlator 
|({S(x^), S(?/^)}), which is defined to only involve measurements near the light cone. In 
the second case, there can be a contribution but only if the measurement points x and 
y are very close to x' and so to each other. That's because y is connected to {x',u) by 
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FIG. 16: Regions of x and y that contribute to the current correlator 2({j'^(^))i'^(y)}) when the 
bulk points, which lie near the light cone, have 4-positions indicated by the black dots. The two 
diagrams correspond to the cases of diagrams (a-1) and (a-2) of fig. 13. 

X y 

I 

>\ 

Xi Xi 

FIG. 17: Labeling of points for discussing diagram (c-1) of fig. 13. 

a low-momentum and so must lie in the forward-time diffusive region relative to x', as 
depicted by the magenta region in fig. 18b, while x is connected by a low-momentum Q,.,. and 
so may lie in either the forward-time or backward-time regions, as depicted by the magenta 
and blue regions of fig. 18b. If both are to be measured close to light cone to determine 
then all the points are close to each other, and so diagrams (c-1) and (c-2) 
only contribute to very localized correlations. Our goal, however, is to understand whether 
^{{T,{x^),T,{y^)}) has support at relatively large separations. To answer that question, we 
may therefore ignore (c-1) and (c-2). 

The low- momentum graviton propagator has a sound pole.^^ Though not necessary for 
the above argument, we note in passing that for a graviton internal propagator in (c-1), 
the picture in fig. 18a therefore plausibly has some relation to the type of physics discussed 
earlier in fig. 7. 

The upshot of this section is that we may compute the refined correlator 



See, for example, ref. [28], which derives sound poles in the bulk-to-boundary graviton propagator and then 
in 2-point boundary correlators. We will discuss later the relation of bulk-to-boundary and bulk-to-bulk 
propagators at finite temperature. 
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(a) 



FIG. 18: Restrictions on where the measurement points x and y can be for diagram (c-1) in the 
cases where the bulk point x' of fig. 17 is (a) far from or (b) near the x^ ~ x^ hght cone. In the 
first cannot be below the lightcone because of causality constraints. 

|({S(2;3),E(?/3)}), of section IID2 with only dia gram (a-1). But the causality prop- 
erties of diagram (a-1) are such that the reasons for carefully defining the refined cor- 
relator no longer apply. So, restricting attention to (a-1), there is no reason not to 
compute the original, simpler correlator i({G(x), 0(?/)}), or its time-integrated version 
|({S0(a;°=oo, a:;), S0(2/°=cxD, t/)}). That is, we will compute (2.35) with the "relevant dia- 
grams" being (a-1). 

IV. THE BULK-TO-BULK CORRELATOR G^^) 

In order to evaluate diagram (a-1) of fig. 13, we will need the bulk-to-bulk correlator 
associated with the internal line of that diagram. (The superscript is a reminder that this 
is the 5-dimensional bulk correlator and not the 4-dimensional boundary correlator.) When 
evaluating correlators of the time-integrated charge deposition Se(a;°=oo, a;), we will turn 
out to need only the case where the two bulk points in diagram (a-1) are at very large times. 
At very large times, the 5-dimensional excitation created by the source has fallen in the 5th 
dimension, so that it is localized very near the horizon. In consequence, we will really only 
need the bulk-to-bulk correlator Gir between near- horizon points. 

In section IV A below, we give a general discussion of how bulk-to-bulk propagators such 

(5) 

as Grr can be computed in terms of more familiar bulk-to-boundary propagators Q. We 
then specialize in section IV B to the case where both bulk points lie very near the horizon. 

A. General case 

For diagram (a-1) of fig. 13, we are interested in bulk-to-bulk propagators for a 
transversely-polarized bulk gauge field A±_ that is dual to the transversely-polarized R- 
charge current operator j^. The generalization to other cases is relatively straightforward. 
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In particular, in this section we will simultaneously treat the case of a bulk scalar field, 
which does not require adding any additional indices or other complications of notation and 
which is relevant to our discussion of ejet 3> 1 in Appendix A. 



1. Bulk-to-bulk G 



(5) 
R 



We start by looking at the retarded bulk-to-bulk propagator G^^ from one bulk point 
Xi — {xi, ui) to another X2 — {x2, U2)- Note that xi and X2 label 4-positions of bulk points 
in this discussion and are not the source points we have labeled xi and X2 previously. Let q 
be the Fourier conjugate of X2—X1. We will verify in a moment that can be written in 
terms of bulk-to-boundary propagators as 



G^i\u,;u2,q) = G^^{u,;u2,q)^ 



G 



WKAiq) 



Gr{u>; q) [Gr{u<; q) - Qk{u<] q)] , (4.1) 



where W^ra(q') is the Wronskian (to be made explicit in a moment), = min(-ui, ■U2), and 
M> = max(Mi, M2). The overall normalization Ag comes from the normalization of the kinetic 
term in the 5-dimensional Lagrangian. For gauge fields, we take >Ckin = ~i(fi'sG-^)^-^^"'"-^/j 
and correspondingly 

^fG = qIgR. (4-2) 

where R is the AdS radius and (jfsG = 47r/A^c [29, 30]. For scalars, >Ckin — 
-{9lGR''r\di<l>)9'\d.j(t>) and Mg = qIgRV^- 

Let's now check that the propagator (4.1) satisfies the equation of motion with a source 
term at {xi,ui): 



1 



where 



A 



g^^^ transverse gauge field Ax_\ 
1, scalar field 0, 

and m = in the gauge field case. In 4-momentum space, this is 



(4.3) 



(4.4) 



1 

'M'g 



d5 {V^ /'A d,) - (Xq^g^'q, + m') G^^ (u^; u, q) ^5{u-u,). (4.5) 



From the fact that the bulk-to-boundary propagators satisfy the homogeneous equation 
1 



ATg 



d5 /'A d,) - ^g {Xq^.g^'^q, + m'')\g{u- q) = 0, 
one can verify that (4.1) satisfies (4.5) provided VFra(?) is the Wronskian 

w^RA = v^^''a(6;J5^a), 



(4.6) 



(4.7) 



which is independent of u. Here, d is defined by fdg = fdg — {df)g. 

For the bulk-to-boundary propagator Q{u\ q), our convention is that q is the momentum 
conjugate to the boundary point. For the bulk-to-bulk correlator G^\ui]U2,q), the con- 
vention is that q is associated with the bulk point at U2, which we take to be the later-time 
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point in defining the retarded correlator. Eq. (4.1) gives a retarded bulk-to-bulk correlator 
because it is analytic in the upper-half g° plane since (i) ^r('u; q) has this property and (ii) 
the poles in ^a(m<; q) cancel the corresponding poles in the denominator Wra- 

Finally, bulk-to-bulk propagators should satisfy the correct boundary condition at the 
boundary, which is that the boundary is fixed and so fluctuation fields (such as (f) and A±) 
should vanish there in this context. That is, G^^\ui;u2,q) should vanish for ui or U2 zero. 
This is easy to see for (4.1) in the case of A± or massless scalar fields, since the bulk-to- 
boundary propagators are then normalized to 1 at the boundary, so that ^a(0; ?) — ^r(0; q) = 
0. Eq. (4.1) also satisfies the boundary condition in the massive scalar case, but there one 
must take care to regulate the boundary. 

2. Bulk-to-bulk 

The advanced bulk-to-bulk propagator is similarly 

Ofiui, U2, q) = G^aH^i; «2, q) = ^^^(^q^ ^a(m>; q) [^r(m<; q) - ^a(m<; q)] ■ (4.8) 

3. Bulk-to-bulk G^? 

Applying (3.3) then produces 

Gifiu,; U2, q) = cth(i^a;) [G'i\u; q) - Gf(u; q)] 

= ctha/3u) — ^ [gK{ui- q) - GAiui; q)] [gK{u2; q) - Qa{u2; q)] . (4.9) 



B. Near- horizon case 

Here and throughout the rest of this paper, we will work in units where 2'kT—I and 
specialize to the choice of 5th dimension coordinate u with metric (2.1). 

1. Near-horizon bulk-to-boundary propagator 

Near the horizon, the general solution to the homogeneous equation of motion (4.6) with 
the boundary conditions appropriate to Qyi can be expanded in powers of \—u as 

g^{u- q) = h{q) (1 - u)~'^/''[l + 0(1 - u)]. (4.10) 

We will not need the exact solution away from the horizon. All of its complexity is 
summarized here in the factor h{q), which is determined by the normalization condition for 



^'^ As an example, for the case of A^, the exact solution is given by replacing [1 + 0{1 — u)] in (4.10) by 
[i(l + u)]"/^Hl(2,g2- 1^2^ izij^. l^uj, l-iw,0; 1-u), where HI is aHeun function as defined 

in ref. [31]. 
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^R, on the boundary (e.g. ^r(0; g) = 1 for A± and massless scalar fields). The answer to 
our question of whether jet charge deposition is correlated over relatively large distances 
will turn out not to depend on the details of b{q) — all that will matter are some simple 
analytic properties of b{q). Since ^r, is a retarded propagator, b{q) must be analytic in the 
upper-half complex frequency plane. Because [^R(xi; —q)]* — Gr{u; q) for real q, (4.10) also 
gives the important relation 

b{-q)=b*{q) for real g. (4-11) 

It will be convenient to change variables from u to 

T~ -iln(l-'w). (4.12) 

The ~ indicates that the details of the definition will not matter with regard to sub-leading 
corrections as w — )■ 1. Physically, r represents the time (as measured by an asymptotic 
observer) for the a wave in the bulk to propagate from near the boundary to the near- 
horizon position u ~ 1, and details such as initial conditions or phase vs. group velocity 
only affect sub- logarithmic corrections to (4.12). In terms of r, (4.10) becomes 

gf{u; q) = b{q) e"^^ [l + ©(e-^^j . (4.13) 



2. Near-horizon hulk-to-hulk 

We now turn to as given by (4.1). For simplicity of presentation, we will specialize 
here to the case of the Ai_ propagator and quote the scalar case at the end. We first need 
the Wronskian (4.7). Since the Wronskian is m- independent, we can evaluate (4.7) at any 

and it is most convenient to evaluate it in the horizon limit u ^ 1. Using (4.10) and its 
conjugate 

Qk{u; q) = GKiu; -q) = b{-q) (1 - u^'^ [l + 0(1 - u)] , (4.14) 

together with 

^ g^'g^^ ^lRf^R{l-u), (4.15) 

we get 

WnA = -tR^b{q)b{-q), (4.16) 

which for real q is Wra = —iRuj\b{q)\'^. Combining (4.16) with (4.1) and (4.13), the near- 
horizon limit of the retarded bulk-to-bulk propagator is 

2 

u 

— ^SG |^_gi%)giw(T2+Ti) _|_ gja)|-n2-Ti|j ^ (4-17) 
OJ 



where 
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T = -r ( \ l^oundary 



T2-1A 



X = oo horizon 



FIG. 19: Schematic representation of the two contributions to the near- horizon bulk-to-bulk re- 
tarded propagator (4.17). 



For real q, the phase 6{q) is real. At poles of in the lower-half complex frequency plane, 
6 becomes —ioo. The right-hand side of (4.17) is finite as oj — because^^ 

6{u, q) ^ linearly as u 0. (4.19) 

For scalar fields, the result is a factor of 2 larger than (4.17). 

To understand the physical interpretation of the two terms on the right-hand side of 
(4.17), consider the metric (2.1) in the near-horizon limit, 

ds^ ~ ^ [-f dt" + dx^ + f rfr^] (4.20) 

with / ^ 2(1 — m) ^ 26^^"^. Note that points with different x are far apart compared to 
points with different t or r. Because of this, the x motion decouples. More concretely, the 
equation of motion (4.5) satisfied by the A_[_ propagator is 

{dl + uj') Gt\n- r, q) ^ 2gl^ 6{r - n) (4.21) 

in the near-horizon limit for fixed q. This looks similar to a flat-space problem in one space- 
time dimension, where the Green function would be proportional to (ia;)~^e*'^'^'^'. Compare 
to (4.17). The second term on the right-hand side of (4.17) corresponds to a signal that 
travels directly from ri to T2, corresponding to a distance \t2 — Ti| in r. The first term, 
in contrast, corresponds to a signal that travels from Ti to the boundary r=0 and then 
reflects back to T2, traveling a total r-distance proportional to T2 + ti. 6{q) represents the 
phase shift that the wave picks up traveling through the region far from the horizon, where 
the dynamics becomes more complicated than the simple asymptotic behavior (4.21). The 
relative minus sign between the two terms is because the wave flips when it reflects from 
the boundary due to the boundary condition that the Green function vanish at the horizon. 
The two contributions to the right-hand side of (4.17) are depicted schematically in fig. 19. 



Presumably there is some simple way to understand (4.19), but we just checked it numerically. 
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(5) 

3. Near-horizon bulk-to-bulk Grr 

Applying the same results to the formula (4.9) for Gr? gives, in units where 27rT = 1, 

CO 

, ctHnuj) I ^ ^ ^(^^^ ^ cos(a;(T2 - n))] . (4.22) 



cos 



Note that the w— j-O limit is finite because of (4.19). 

The second term in (4.22) is independent of q and so is a 5-function d^^^Xi — X2) in 
3-position space. This locality is a consequence of the fact that x motion is suppressed near 
the horizon. In contrast, the first term docs depend on q through the phase 5{q) and so can 
produce non-local correlations in x. This non-locality arises because the phase 5 appears 
in contributions like the dashed line in fig. 19, corresponding to propagation away from 
the horizon, so that 3-space dynamics are then nontrivial, bouncing off of the boundary, 
and then propagating back to near the horizon. We will see that only the second term in 
(4.22) contributes significantly to charge deposition correlations for high-energy jets, and the 
locality of that term will be responsible for the locality of the charge deposition correlation. 

In Appendix B, we show how to relate our result (4.22) to expressions derived for bulk- 
to-bulk Grr by Caron-Huot, Chesler, and Teaney [10] in the context of fluctuations of a 
classical string. 

V. CALCULATING THE CORRELATOR 
A. Starting formula 

We are now finally ready to calculate the contribution of diagram (a-1) of fig. 13 to our 
correlator. In our application, only the standard Yang-Mills 3-point vertices contribute to 
this diagram: the contributions from the Chern-Simons term in the SUGRA Lagrangian 
vanish because the SU(4) group factor d^^^ = 2 tr [{T^*", T^}T'^] that would be associated 
with a Chern-Simons 3-point vertex vanishes. 

To save space, we will henceforth abbreviate {\{3^^'^''{x), j^^'^''{y)}) as {^{j^J"}). Plug- 
ging the source (2.5b) into the relationship (3.11) between and the 4-point equi- 
librium correlator Gaarr, and expressing Gaarr in momentum space, gives 

yV'2 /• 

A(|{j^ f}) = ^ ^ / e^ep G^^f'^'^f^^'^iQi, Q2; Q, Q') A^Qi - k) KHQ^ + k) 

^ JqQ'QiQ2 

X e^e-^e^^'-^(27r)^5(^)(g+g'+gi+g2). (5.1) 
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(a-1) 

FIG. 20: Naming conventions for momenta assignments to diagram (a-1). 
Labeling momenta as in fig. 20, diagram (a-1) plus its permutations gives 

' du,{^gg^^gniQt{u,-Q,)g^^{u,-Q) 



i{-p-Qi)p 

X 



p=Qi+Q=~{Q2+Q') 

+ {Q,fi^Q\iy}. (5.2) 

The group factors are /"'^^/^''"^ = —2. In Appendix C, we show that these formulas give a 
simple and easily predictable result for A{^{j^, j'^}) in the T=0 case. But here we focus on fi- 
nite temperature and the contribution of diagram (a-1) to ^{{f^Q{x^=oo, x), T,Q{y'^=oo, y)}), 
which we will abbreviate as |({Se, Se)} in the following. 

From the definitions (2.16b) and (2.15) of T,q and O, replacing and by Ee in (5.1) 
is equivalent to replacing 

^.^'^(^i; Q) ^ '^''qo^^' Q) (5.3) 

and similarly for Qf-" {u2] Q') in (5.2). Note that the R-charge diffusion constant (2.14) is 
D=l in our units 27rT=l. As reviewed earlier, our measurements do not need to resolve 
scales as small as 1/T in order to study jet stopping, and so we may make small-momentum 
approximations for the momenta Q and Q' conjugate to the measurement points x and y. 
In previous work [3], we discussed how the corresponding low- momentum bulk-to-boundary 
propagators are 

k) - - -^(1 - uY-^-/^ (5.4) 

ZLO — k'^ too — k"^ 

k) - -- ^ + -^(1 - u)-'\ (5.5) 
luj — k^ luj — k^ 

We also showed that for the purposes of computing {T,q{x^=oo, x)), only near-horizon values 



30 



of u were important and that in that hmit one could replace^^ 

(5.6a) 

and ignore altogether in the calculation: 

~ 0. (5.6b) 

The same will happen in computing the correlation of Ee- Using (5.6) in (5.3), we see that 
replacing the j's by Ee's in (5.1) amounts to replacing 

g^^{ui-Q)^5l and g^"" {u2,Q') ^ 5^ (5.7) 

in (5.2). Putting everything together and using 



gives 



(|{Ee,S0})(a-i) - 7^ / A*{x,ui)d^oiGf^^{x,ui]y,U2)dyoA{y,U2) 

^9SG J ""l/l ""2/2 



J\f^ r dui du2 
+ {x^ y}, (5.9) 



where 



A{x, u)= j gl{u; q) Mq - k) e"^-^ (5.10) 

Jq 

describes the bulk vector field created by the source and 

A*{x, u)= I g^{u- q) Kl{q - k) e-"^ \ (5.11) 

Jq 

For comparison, the corresponding formula for (Se) is [3] 

(Ee)~4^ / — A\x,u)id,.A{x,u). (5.12) 

At this point, it is useful to think again in momentum space, where the operators 5^.0 
and dyf) in (5.9) become the —i{—p — Qi)o and —i{Q2 — p)o factors in (5.2) respectively. As 
mentioned earlier, our resolution requirements permit approximating Q and Q' as small, in 
which case these factors may be approximated by 2i(Qi)o and — 2i(Q2)o- Because of the 
source factors A^, in (5.1), our calculation only gets contributions from {Qi)o — —E and 
{Q2)o — E. The upshot is that the two operators do in (5.9) may each be replaced by —2iE: 

(i{SG,Se})(a-i)^-2£;'^Re / ^ ^ A*ix,m)iGil{{x,m;y,U2) A{y,U2). (5.13) 

9sG J ■"1/1 ■"2/2 



See specifically the discussion in sections IV.A.l and IV.G of ref. [3]. 
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B. Factorizing G^^^ 



Let us return to momentum space for G 



(5). 
rr ■ 



(i{Se,M)(a-i)^ 



2E'^Re [ [ ^^A*{x,u,)e'P-'tG'^{{u,-U2,-p)e-'PyA{y,U2). (5.14) 

fl'SG JpJ '"'l/l '"'2/2 



Using the general formula (4.9) for Or? and the result (4.16) for VFra, the ui and U2 integrals 
factorize as 

A/"! „ f cth(7r/) 

where 



;HS0, Se})(a-i) - ^ Re J*ip, x) J {v. v) , (5.15) 



J{p,x)^ J K(«;p)-^^(«;p)]e-^-^(a;,«) 

= ^ A^g - ^) e^(^-^)- J (i.; g) [^^(i.;p) - gj:(u;p)] . (5.16) 

C. The u integral 

We now evaluate the u integral in (5.16), 

= y gliu;q) [g^iu;p) - ^^(«;p)] . (5.17) 

Consider the near-horizon (large r) approximation to the integrand. Using (4.13) for 
and using du/uf ~ o?t/^ 

/•oo _ - 

Iqp ~ b{q) J dr e^«°^ [b{p)e'P°^ - b*(p)e-'P°^^ . (5.18) 

Recall that g° is really q'^ + ie because of the retarded prescription on q. The integral is then 

b{p) h*ip) 



(5.19) 



_gO+p° + ie q^-p^ + i€_ 

Why is it okay to make the near-horizon approximation for the bulk-to-bulk propagator? 
We are interested in evaluating the correlator (5.14) of Se's at x^,y^ — ?> oo. In this limit, 
A{x,Ui) and A{y,U2) are exponentially suppressed everywhere except very close to the 
horizon (r > 1) [3].^^ 



^"^ We will consider only real values of p (but complex values of q) in what follows, and so b{—p) and b*{p) 
are interchangeable as in (4.11). 

We can see this in detail from the exponential tail derived in our previous work [3] for large X+ = x+— t(u): 
A{x,u) ~ ie'm=o-+r{u)] j^gg A^^\q^^^ ; X- +t{u)) e" . 

Large a;° for fixed corresponds to large —x~ and x^ . From the x~ + t argument of A^, we therefore 
get exponential suppression unless r ~ —x~ is also large. In this case, X+ ~ 2x^. 
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Note from (4.13) that corrections to the near-horizon Hmit formulas for our propagators 
are suppressed by e"^"^, and so these corrections vanish in the large-time limit of interest. 

D. The g° integral 

Putting the result (5.19) for the u integral back into eq. (5.16) for J gives 

b{q)~A,{q-k) J,.,. -ib*{p) h{q)K,{q-k) ' . , , • (5-20) 
q -r p -r It jg q p -r it 

Now do the q^ integral by closing the integral in the lower half-plane/^ and look at the 
contribution from picking up the explicit poles above: 

J{p, x) ~ e^^f°^°6(p) f h{-p^, q) Al{-P^ -E,q-k) e^(<'-P>=" 

Jq 

- b*{p) I h{p\ q) Al{p' -E,q-k) e'^'^-^^^ . (5.21) 

Jq 

Before we use this expression, wc should consider what other singularities contribute to the 
integral. In particular, h{q) has poles in the lower-half q^ plane, corresponding to quasi- 
normal modes. When evaluated at these poles, the e~'^°^° term in the integrand will give a 
suppression factor of order 

e^"^'"''p°i«. (5.22) 

(Keep in mind that Im ^'p^^j^ is negative.) But we are interested specifically in the x° — > 00 

limit for evaluating Eej and so these contributions will vanish. Because Kl only has support 
for g_ ~ where the poles in q^ have imaginary parts of 0{E~^^^) [3], the relevant poles 
in q^ = q_ — q^ have imaginary parts of —0{E^^^'^). 

Because of the support of the envelope factor A^, in the first term of (5.21), the e^'*'°^° 
phase factor in that term has p° ~ and so is highly oscillatory. Remember that we are only 
trying to resolve distance and time scales of the observables on scales large compared to 1 /T. 
If we smear out our observables "Eeix) and Se(y) over such scales, then the contribution of 
such a highly-oscillating phase will be smeared away.^° As a result, we may drop the first 



Technically, some care should be taken here because of the behavior of Aj, in the complex frequency plane. 
Really, one should take care to route the closing piece of the contour through regions of the complex plane 
where the integrand is exponentially suppressed. See Appendix E of ref. [3] for related discussion (though 

there the contour is different and serves a different purpose). 
More specifically, suppose we replace T,e{x) by 

[Ee(a;)]snieared = d {Ax) T,0{X + Ax) -2 

J ''^smear 

as discussed in sec. II B of ref. [3], where the smearing distance ^smoar is chosen large compared to 
microscopic scales such as 1/E and 1/T but small compared to the scales we're interested in resolving, 
such as stopping distances. If one applies this procedure to something that behaves like exp{i2Ex^), then 
one obtains an exponentially small result. 
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term of (5.21) and simply write 

J(p, x) ~ -b*{p) I h{p\ q) Al{p' -E,q-k) e'^'^-^^^ . 

Jq 



(5.23) 



E. Assembling the pieces 

Now use the expression (5.23) for J' in the correlator (5.15) of Ee's: 
(H^e,Se})(a-i) ^2£;^:^Re / S^M^e^P-M f ^-i'.-^^-W v 

9SG Jp P Jqq' 

X 6*(p°, q) 6(p°, q') Al{p'-E, q-k) Mp'-E, q'-k). (5.24) 
Now note that the p integral just gives S^^\x — y): 

(i{Se,Ee}>,.-i,-2£^^«'«(a^-»)Ro/ "^^^ f 

9sG Jp" P Jqq 



9SG 

X 6*(p°, q) h{p\ q') Alip'-E, q-k) Al{p'-E, q'-k). (5.25) 

The right-hand side vanishes for x ^ y. Recall that throughout we have made approxima- 
tions that blur our resolution of x and y on the scale 1/T. 

And so we have answered the question we set out to answer: For \x — y \ large compared 
to 1/T, the correlation is exponentially small (coming from all the exponentially-small cor- 
rections that we dropped throughout). 
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Appendix A: Generalization to jets with large R charge 

In this appendix, we discuss how the argument in the main text is modified if one uses 
a source operator with R charge larger than one. For simplicity, we will consider scalar 
operators. An example is tr(X^), where X is any one of the three complex scalar fields 
in Af=4: SYM. This operator has conformal dimension A. The R charge J of tr(X^) is A 
under a U(l) subgroup of the SU(4) R charge symmetry, which is the subgroup that we will 
choose for our measurement operators j°(x) and j^{y). [That is, it is the U(l) subgroup 
that corresponds to t^/2 in the main text.] 

In what follows, we more generally consider taking our source operator to be any scalar 
BPS operator of conformal dimension A and R charge J under a U(l) R symmetry subgroup. 
In the 5-dimensional bulk these correspond to scalar fluctuations of mass m such that 

A = f + ^(f)2 + m2i?2, (Al) 
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where d=4 is the number of boundary space-time dimensions. The important property of 
the bulk-to-boundary scalar propagator is that it takes the form 

gn(u, q) = b(q) (1 - u)-'^/^[l + 0(1 - u)] = b(q) e'^^[l + ©(e"^^)] (A2) 

near the horizon, similar to (4.10) and (4.13) but with a different function b{q) than in the 
transverse gauge boson case. Near the boundary u—0 it behaves like 

gniu,q)o,u^'-''y'[l + Oiu)]. (A3) 

The divergence of the bulk-to-boundary propagator at the horizon for d > A (related to the 
need to renormalize such operators) will not have any effect on our calculation, which will 

be controlled by the near-horizon behavior. 

In the vicinity of the horizon, the bulk-to-bulk retarded and symmetrized propagators 
take the form 

2 

iCiJ^ (n ; T2 , g) = 2 ^e^'^^ [-e^-^^^) e^'^^< + e-^'^^< ] (A4) 

00 

iGif{n;T2, q) = igla ^^^^^^^ J_ cos (^cu (t2 + n) + 5(g)) + cos(^a;(r2 - ri))] , (A5) 

where, as before, e*'^'^'^^ = h{q)/h{—q). These formulas differ from the transverse boson case 
(4.17) and (4.22) by a factor of 2 normalization (and 5{q) is a different function of q). 

For a scalar field with R charge J, the same real-time Witten diagram (a-1) as in the 
main text gives 



I (a-1) 

X 



X iG'[^)(Mi;K2,-p) 

P=Qi+Q=-{Q2+Q') 

+ {g,/xf^g',i/}. (A6) 



The full expression for the scalar bulk-to-boundary propagator can be written as 
where 

h{u) = Hl(2, q^ + ^- iA_uj + ii±^ - ^ + ^ + 1 - iw, A_ - 1; 1 - u), 



-|- rn^E? = d — A, and HI is the Heun function [31]. (Compare to footnote 14.) Here we 
have naively normalized the bulk-to-boundary propagator so that the coefficient of the small-u behavior 

(4u)^-/^ = (47i)('^^'^)/^ = z'^^^ is precisely 1 instead of some more complicated normalization involving 
regularization of the boundary to some small, finite vahic zq of z = \/4m. As argued by ref. [30] , the naive 
normalization of the bulk-to-boundary propagator can be used when computing n-point functions with 
n > 2. 
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Substituting the metric factors and the various propagators leads to 



/I ^ TV A/lJ f dUi dU2 .^fsW ^ 

(2|Ee,E0})(a-i) ~ — ^ / A {x,ui)d^oiG\,>{x,ui]y,U2)dyoA{y,U2) 



where now 



/SG "-1^1 '^2/2 



(A7) 



(A8) 



The same considerations we have previously made apply: namely we can approximate the 
action of and dyO by —2iE factors, substitute the factorized form of the symmetrized 
propagator, and evaluate the subsequent u integrals. 



"■QP 



du 
IF] 



(A9) 



Iqp above differs from (5.17) by a factor of 1/-U, but in the near-horizon approximation ~ 1 
and du/{v?f) ~ dr it becomes formally identical to (5.19): 



J-qp — 



ib{q) 



h{p) 



h*{p) 



(AlO) 



The final result is then the same as (5.25) but with the substitution M\ — )■ M\J'^/A in 
the overall normalization. The conclusion is the same as for the case analyzed in the main 
text: for late times {x^ , ^ 00) the integrated charge deposition correlator is exponentially 
suppressed for \x — y \ ^ 1/T. 



Appendix B: Relation of our near-horizon to ref. [10] 

Caron-Huot, Chesler, and Teaney [10] studied correlators of fluctuations of the position 
,x(t, r) of a classical string, where r was their coordinate for the 5th dimension, running from 
r=l at the horizon to r=oo at the boundary. One of their results was 

G'rr(ii,ri;t2,r2) = j dt[ dt'^ [-G,^{ti,n;t[,r^)] [-G,^{t2,r2;t'^,ri,)]G^{t[,Q (Bl) 

for ri,r2 > rh, where rh = 1 + e with e very small defined a "stretched" horizon and where 
G^^{t'i, ^2) was a type of correlator on that stretched horizon. For our purposes here, we may 
treat e as infinitesimal. They found that 

G^ = --d^,d^. In |1 - e-2-^(''i-''2)|, (B2) 

TT 

where v is Eddington-Finkelstein time 

1 f dr ,^ , 

'^*^^J7^ (B3) 

evaluated in (B2) at ri = r2 = rh. So (B2) can be rewritten as 

Gl = -"^dtA, In |1 - e-^'^^^*!-*^)!. (B4) 
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Their definition of Grr as |({0(1), 0(2)}) for an operator O is i/2 times ours, but we will 
not worry about overall normalization factors in this discussion since the final formulas for 
Grr also depend how one normalizes the operators of interest to a particular problem. 

We will now see that these same equations describe formulas in our paper (up to overall 
normalization) provided we generalize (B2) to fields with x dependence, 

G^^\xi,n]X2,r2)= / dx'^dx'^ [-G'(^)(a;i,ri;a;;,rh)] [-G'(^na^2,r2;4,^h)]G!;(a;;,4), (B5) 



and keep (B4) the same. 

Fourier transforming (B4) gives 

G),{q)o^ujci\i{\^uj). (B6) 

The Fourier transform of (B5) is 

= [G^i\n.;n,q)rG^i\r^;r2,q)G'M, (B7) 

given our convention in this paper that Gr means Gar- 

We can now proceed generally or specialize to the case where ri and r2 are near the 
horizon (though not as near as rh). In the latter case, simply plugging in our near-horizon 

formula (4.17) for gJJ^ and (B6) into (B7) gives 

which reproduces our near-horizon result (4.22) for GiiK In the more general case, use (4.1) 
instead of (4.17) to obtain 

Gif oc u cth(i/3a;) {^^^ [G^im; q) - Q^UV, q)] }* {^^^[Gniu^; q) - Q^{u2; q)] ] 

This reproduces (4.9) with the aid of (4.10) for ^R(wh; ?) and (4.16). 
Appendix C: T—0 analysis of j"}) 

It's useful to see how the formalism used in this paper works in the zero-temperature 
case. Among other things, it provides a useful double check of the normalization of some of 
our basic equations. 

There is no thermalization and "charge deposition" in vacuum, and so wc will study 
^{{j^{x),j'^{y)}) instead of |({Se(^)i ^eiv)})- At zero temperature, a jet excitation simply 
propagates forever. The excitation's total current should not change with time, and so it is 
natural to expect that 

l^{\{jj>'{x\x),jj{y\y)])^^^^[jj^{x\x))^^(^jj\y\ 

^el'^e{x')e{y^) (CI) 
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at zero temperature. The goal of this appendix is to show that the methods used in this 
paper indeed reproduce this expected resuh. Here, ejet = 1 as in the main text. Above, the 
A indicates we subtract (vacuum) fluctuations as in (3.6), analogous to (2.31). 
In the following, we will write the zero-temperature metric as 



4 



2 



-(-df + dx"^) + ^dU 
u w' 



(C2) 



where u corresponds to our earlier ■u/(27rT)^ and runs from zero to infinity. 

Start from (5.1) and (5.2). As discussed in ref. [3], the low-wavenumber approximation 
for the observable corresponds to 

gf^{u;Q)c:.6ii (C3) 
at zero temperature. In this approximation, (3.14) gives 

g^'^{u;Q)c^O, (C4) 

and so we can ignore every diagram in fig. 13 except (a-1). Combining (5.1) and (5.2) using 
(C3) gives 



A(Hi^n)-f^ ——A*(x,u,)d,,iGll[{x,u,-,y,U2)dy.A(y,U2) 

+ {x^ y}. (C5) 



Note that the zero-temperature formula for happens to have exactly the same 

form as the finite-temperature formula (5.9) for |({Se,Se}) (noting that /— >1 at zero 
temperature). Following the same line of argument that led to (5.14) then gives the analogous 
zero-temperature formula 



A(i{/,n) 

-2^'^r^Re / I :^:^A*{x,ui)e'P-''tG^l>^{uuU2,-p)e-'P-yA{y,U2). (C6) 

It's helpful at this point to have the explicit formulas for the transverse-polarized bulk-to- 
boundary propagators in 4-momentum space, which are 



giu; q) = yW^i( V4^), (C7) 

where Kn is the modified Bessel function of the second kind, in this formula means q^ + ie 
for the retarded propagator and g° — ie for the advanced propagator. The Wronskian (4.7), 
which is most easily evaluated in the it ^ limit, is then 

iyRA(?)=M?'sign(5'^). (C8) 

Note also that the zero-temperature (/3 — )■ oo) limit of the fiuctuation-dissipation relation 
(3.3) is 

iGrriq) = sign(?°) [iCM - iGA{q)] ■ (C9) 
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Comparing (C8) to the finite-temperature formula (4.16), the Tipshot is that the analysis 
of at zero temperature goes through just as the finite-temperature derivation of 

(5.15) for (^{Se, but with the replacement 
to give 

A(i{j^r})^-2Pr:^Re / Ji^r{p,x)J{p,y), (Cll) 

with J" defined as before in (5.16). Using (C7), the desired u-integral (5.17) of J is given 
explicitly by 

POO 

Iqp = Am ei-p") ^/-p^q^ / duK^i ^W) Ji ( ^ -Aup^) 

Jo 



^2m^--J{-p') (C12) 



p^ — q 

with q'^ interpreted as 5° -|- ie. Substitution into (5.16) gives 

Jip, x) = -2mp^ 0{-p') / Ai(? - k) (C13) 

Jq r - 

as the zero-temperature analog of (5.20), but note that we have not made anything analogous 
to the "near-horizon" approximation in the zero-temperature analysis. 

Since our sanity check (CI) of results and normalizations in this appendix involves inte- 
grating over 3-position, it's convenient at this point to integrate (C13) over x to get 

Jj^ ip'e{-p^) J dq^ Uq^-E^p\p^-E). {CIA) 

Following similar approximations as from (5.20) to (5.23), pick up only the explicit poles in 
(C14) and then throw away highly-oscillatory terms in the result. This yields 



J~^^(V)^(^°)Al(p-^) 



~ 47rp+ e{-p+) e(x°) Kl{p - k). (C15) 
Using this expression in (Cll) gives 

A(H/.j"Jyn)^-327rfc''r^(x°)^(y°):^ / ^(-P+)^ Al(p-^) 

fl'SG Jp P 

^Q^e{x')e{y% (CI6) 

where 

q-SttE^ / e{-q^)\q^\\KL{q)\'. (C17) 

fi'SG Jq 
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This formula for Q is equivalent to that found in ref. [3] for the average charge created by 
our source operator. Dividing both sides of (C16) by Q and invoking (2.19) finally gives us 
(CI), as expected. 
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